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The central theme of complex systems research is understanding the emergent macroscopic prop- 
erties of a system from the interplay of its microscopic constituents. Here, we ask what conditions 
a complex network of microscopic dynamical units has to meet to permit stationary macroscopic 
dynamics, such as stable equilibria or phase-locked states. We present an analytical approach which 
is based on a graphical notation that allows rewriting Jacobi's signature criterion in an interpretable 
form. The derived conditions pertain to topological structures on all scales, ranging from individual 
nodes to the interaction network as a whole. Our approach can be applied to many systems of 
symmetrically coupled units. For the purpose of illustration, we consider the example of synchro- 
nization, specifically the (heterogeneous) Kuramoto model and an adaptive variant. Moreover, we 
discuss how the graphical notation can be employed to study isospectrality in Hermitian matrices. 
The results complete and extend the previous analysis of Do et al. [Phys. Rev. Lett. 108, 194102 
(2012)]. 



Characterizing the behavior of complex systems and 
discovering the critical boundaries in parameter space at 
which qualitative changes occur are of central interest 
in statistical physics. Many studies have demonstrated 
that, in specific models, this can be achieved by the in- 
struments of nonlinear dynamics [D-Ql- However, apply- 
ing these tools to large heterogeneous systems, such as 
complex networks of dynamical units poses significant 
challenges (B-Q- I n particular, there is no simple way 
to capture the structural aspects of such systems, which 
are known to crucially influence the emergent dynamics 

The interplay of dynamics and structure in networks 
has extensively been studied using the example of syn- 
chronization yl, [ll|, [l2[. Systems of coupled phase os- 
cillators, originally proposed by Kuramoto have 
become a paradigmatic example. The exploration of such 
systems revealed that the synchronizability of a network 
depends on global topological measures such as the clus- 
tering coefficient, the diameter, and the degree or weight 
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but also on local and mesoscalc de- 
This shows that a conceptual under- 
standing of collective phenomena on dynamical networks 
can only be achieved if structural properties on all scales 
are taken into account. 

In a recent paper, we proposed a nonlinear dynamics 
approach, which allows to study the dependence of dy- 
namics on structure in networks of symmetrically coupled 
units [13, [H|. We showed that Jacobi's signature crite- 
rion (JSC) can be used to determine necessary conditions 
an interaction network has to obey in order to support 
a given type of macroscopic behavior. These conditions 
pertain to subgraphs on all scales, from single nodes to 



the entire network. In this paper, we complete and ex- 
tend the analysis of [H|. We present the full deriva- 
tion of the proposed approach and discuss a number of 
new applications including an adaptive Kuramoto model, 
in which the topology of the coupling network coevolves 
with the local dynamics of the phase oscillators J2(| H3] ■ 



I. STABILITY IN NETWORKS OF 
PHASE-OSCILLATORS 

The approach outlined below is applicable to all dy- 
namical systems with Hermitian Jacobian matrices, and 
in particular to all networks of symmetrically-coupled 
one-dimensional dynamical units. However, for the sake 
of illustration, we first consider specifically the example 
of N coupled phase oscillators 
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Here, Xi and w, denote the phase and the intrinsic fre- 



quency of node i, and A £ 
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is the weight matrix 



* ly@mpipks-dresden.mpg.de 



of an undirected, weighted interaction network. Two os- 
cillators i,j are thus connected if Aij = Aji ^ 0. Equa- 
tion (rj} represents the so-called Kuramoto model, that is 
today considered to be a paradigm for the study of syn- 
chronization phenomena in coupled discrete systems 14 1 . 
The model is therefore used as the natural benchmark for 
comparative evaluations of performances of methods and 
tools. 

It is known that the Kuramoto-model can exhibit 
phase-locked states, which correspond to steady states 
of the governing system of equations. The local stabil- 
ity of such states is determined by the spectrum of the 
Jacobian matrix J £ Ht NxN defined by = dxi/dxk- 
As J is symmetric, its spectrum is real. The state under 
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consideration is asymptotically stable if all eigenvalues 
are < 0. 

Below, we analyze the spectrum of J by means of Ja- 
cobi's signature criterion (JSC). The JSC (also known as 
Sylvester criterion) states that a hermitian or symmetric 
matrix J with rank r = N — 1 has r negative eigenvalues 
if and only if all principal minors of order q < r have the 
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Here, the principal minor of order q 
det (Jik), i,k = si, s„. 



Stability analysis by means of JSC is well-known in 
control theory |28f and has been applied to problems of 
different fields. For instance, the JSC has been used in 
[29l | to analyze advection-diffusion equations, and in [30[ 
to determine stability boundaries in thermodynamic sys- 
tems. Moreover, it has been applied to study the dynam- 
ics of offshore platforms [3l| and social systems [32| . 

Two obstacles that are encountered in the application 
of the JSC to large networked systems are the growth of 
(a) the number of terms in each determinant, and (b) the 
number of determinants that have to be checked. Diffi- 
culty (a) can be overcome by means of a graphical nota- 
tion, which allows to derive topological stability criteria. 
This notation is introduced in Sec. [TTJ Difficulty (b) can 
be evaded by considering necessary rather than sufficient 
conditions for stability: While the sufficient condition 
for stability necessitates the evaluation of all principal 
minors, a necessary condition for stability is already ob- 
tained by demanding sgn (D q ) = (— f ) q for some q. 

The necessary stability condition that is found by con- 
sidering a principal minor of given order q depends on 
the ordering of variables, i.e., the ordering of rows and 
columns in the Jacobian. By considering different or- 
derings the number of conditions obtained for a given q 
can therefore be increased [32j. To distinguish minors 
that are based on different orderings of the variables, we 
define S = {si, . . . , s q } as a set of \S\ = q indices and 
D 1515 as the determinant of the submatrix of J, which 
is spanned by the variables x si , . . . , x s . Therewith, the 
conditions for stability read 



sgn (£>| S | )S ) = (—1)1^1, VS, |5| = f,. 



II. GRAPHICAL NOTATION 
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Below, we address difficulty (a) mentioned above, i.e., 
the combinatorial explosion of terms that are needed to 
write out the conditions for increasing \S\. In the com- 
mon notation, more than 700 terms are necessary for 
expressing the minors of order 6. For circumventing this 
problem we employ a graphical notation that captures 
basic intuition and allows for expressing the minors in a 
concise way [24[ . 

The graphical notation relies on a topological reading 
of the minors. We interpret the Jacobian J as the weight 
matrix of an undirected, weighted graph Q. A Jacobian 
element Jij then corresponds to the weight of a link con- 
necting nodes i and j. Moreover, we relate products of 
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FIG. 1. Example for the graphical notation. Shown are the 
minors of the matrix ((Sj in algebraic, and graphical notation, 
and, for each term, the corresponding subgraph of a three- 
node graph Q. Here, as well as in the next figures, filled 
symbols correspond to nodes £ S, open symbols to nodes 

is. 



the Jacobian elements to subgraphs of Q spanned by the 
respective links. For instance JijJjk is interpreted as 
path i — j — k. Similarly, JijJjkJki is interpreted as a 
closed path from i to j to k and back to i. 

We can now express the minors of J as sums over sub- 
graphs of Q. The Leibniz formula for determinants (33j 
implies that (i) a minor D \s\,S is a sum over \S\l elemen- 
tary products Ji 1 j 1 • . . . • J<| S ij| S | ; and that (ii) in each of 
these products every index s; € S occurs exactly twice. 
In the topological reading this translates to the follow- 
ing statements: Because of property (i), each term of a 
minor D^\,s corresponds to a subgraph with \S\ links. 
Because of property (ii) , these subgraphs can be decom- 
posed in cycles of Q: Every index Sj £ S occurs either 
with multiplicity two on a diagonal element of J, or, with 
multiplicity one, on two off-diagonal elements of J. In the 
former case, the respective factor corresponds to a self- 
loop of G, i.e., to a cycle of length n = 1; in the latter 
case, there is a set of factors Jjj, i ^ j, corresponding to 
a closed path of links, i.e., a cycle of length n > 1. 

In summary, the index structure of every term of a mi- 
nor D\g\,s can be expressed by a combination of symbols 
denoting cycles of Q of a given length n. The idea is 
now to supplement the basis of symbols with a summa- 
tion convention; This convention is designed such that all 
algebraic terms that are structurally identical and only 
differ by index permutations can be captured in one sym- 
bolic term, which drastically reduces the complexity of 
the expressions. 

Below, we use the following definitions: The basis of 
symbols is given by x , |, A, □, O, . . . denoting cycles of 
length n = 1,2,3,4,5,.... The summation convention 
stipulates that in a minor D\g\g, every product of sym- 
bols denotes the sum over all non-equivalent possibilities 
to build the depicted subgraph with the nodes G S (cf. 
Fig. 1). With these conventions the first 4 principal mi- 
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nors can be written as 
Di,s = x 

£>2,S = X • X - I 

£>3,s = x- x- x- x- |+2A 
Da.s = x • x • x • x — X • x • I + j 



(3a) 
(3b) 
(3c) 

| + 2 x -A -2n 

(3d) 



More generally, 

D\s\,s = a U combinations of symbols with n = l^l 

(4) 

where symbols with n > 2 appear with a factor of 2 that 
reflects the two possible orientations in which the corre- 
sponding subgraphs can be paced out. Symbols with an 
even (odd) number of links carry a negative (positive) 
sign related to the sign of the respective index permuta- 
tion in the Leibniz formula for determinants [331 ] . 

An example for the graphical notation is presented in 
Fig. [TJ The figure displays the three principal minors of 
the symmetric 3x3 matrix 



Jll J 12 Jl3 
J= | Jl2 J22 J23 

J 13 J23 J33 



(5) 



in algebraic, and graphical notation. Moreover, it dis- 
plays for each term the corresponding subgraph of a 
three-node graph Q . 



III. ZERO ROW SUM 

In many systems, including the standard Kuramoto 
model, fundamental conservation laws impose a zero-row- 
sum condition, such that Ja = — Jij ■ Using this 
relation, we can remove all occurrences of elements Ja 
from the Jacobian and its minors [24| • In the topological 
reading, this substitution changes the graph Q by replac- 
ing a self-loop at a node i by the negative sum over all 
links that connect to i. 

The simplification of the minors due to the zero-row- 
sum condition can be understood using the example of 
the Eqs. ([3]). Replacing the self-loops, the first term of 
every minor D\g\g, x' s l, is (— l)' s times the sum over 
all subgraphs meeting the following conditions: (i) the 
subgraph contains exactly \S\ links, and (ii) it can be 
drawn by starting each (undirected) link at a different 
node in S. By means of elementary combinatorics it can 



be verified that all other terms of D\g\,s cancel exactly 
those subgraphs in x' 5 that contain cycles. This enables 
us to express the minors in another way: Defining 

= y , au acyclic subgraphs of Q meeting (i) and (ii) 

(6) 

we can write 
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FIG. 2. Symbolic calculation of a minor using the zero-row- 
sum condition. Shown is the graph Q, defined by the off- 
diagonal entries of Eq. ([8|. The terms of the minor Di t s can 
be written asx-x-x-x=A + B + C + 2D, — X • X • | = 
-(B+2C), | | = C, 2x- A = -2D and -2a = (cf. Eq. pd])). 
It thus follows that D4.S = $4,3 = A is the sum over all 
acyclic subgraphs of Q meeting conditions (i)-(iv). 
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FIG. 3. Topological equivalents of minors. Shown is the com- 
plete sequence of $\s\ Si \S\ = 1, ■ • • , 5 for the graph Q from 
Fig. El 



We remark that for \S\ = N — I, <£>|s|.s is the sum over 
all spanning trees of Q (cf. Eq. ([5])). Thus, Kirchhoff's 
Theorem [34| . which has previously been used for the 
analysis of dynamical systems (35l . I3H ] , can be regarded 
as a special case of the more general Eq. (0 • 



The simplification of the minors due to the zero-row- 
sum condition as well as the relation between the D\s\,s 
and their topological equivalents $|s|.s can be illustrated 
by means of a simple example. Consider the symmetric 
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zero-row-sum 6x6 matrix 



/(Jll) Jl2 
Jl2 
Jl3 



J 
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(J22) J23 

J23 (J33) J34 

J34 (J44) J45 

J45 (J55) ^56 

J 56 (J 66 )/ 



(8) 



Here, (Ju) abbreviates —Y^jM^ij- ^ n Fig. HJ we calcu- 
late the minor D i} s={i,...,4} in terms of the subgraphs of 
the corresponding graph Q. The calculation illustrates 
the reasoning that lead to the Eqs. (|6]) and ([7]). 

On the one hand, the complete sequence of minors 
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\S\ = 1, . . . , rank(J) can be calculated as 
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where Sk = {1, . . . , fc}, and rank(J) = 5 due to the zero- 
row-sum condition. 

On the other hand, we can use the definition © to con- 
struct the sequence ^isi.s, \S\ = f , . . . ,rank(J), directly 
from the graph Q (cf. Fig. [3]). A comparison of both, the 
algebraic and the topological results, reproduces Eq. (0. 
We note that labeling the nodes in different order would 
have yielded different algebraic as well as topological ex- 
pressions. 



IV. TOPOLOGICAL STABILITY CONDITIONS 

Let us shortly summarize what wc obtained so far. The 
topological reading of determinants maps a Hermitian 
Jacobian J with zero row sum onto a graph Q, whose 
weight matrix is given by the off-diagonal part of J. The 
minors of J can then be interpreted as sums over terms 
associated with subgraphs of Q. Combining the Eqs. @ 
and Q , the algebraic stability constraints on the minors 
of J translate into 



$|5|, s > 0, VS, \S\ = l,...,rank(J). 



(9) 



We emphasize that the graph Q is not an abstract con- 
struction, but combines information about the physical 
interaction topology and the dynamical units. For ex- 
ample, if a graph Q has disconnected components, there 
is a reordering of the variables Xi, such that J is block 
diagonal. This implies that the spectra of different topo- 
logical components of Q decouple and can thus be treated 
independently. 

From Eq. @ we can immediately read off a weak suf- 
ficient condition for stability: Because $|s|,s is a sum 
over products of the J,,, a Jacobian with J™ > Vi,j 



Let C be an unbranched segment of a cycle of Q that consists 
of d links Ci, i 6 {1, . . . , d}. Let further Wi denote the weight of 
Ci, and c x denote the only link in C with negative weight w x . 
Then stability requires that 
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FIG. 4. Stability sets an upper bound on the absolute value 
of a negative link weight. The bounding relation is derived in 
Appendix B. 



is a solution to Eq. © irrespective of the specific struc- 



ture of Q [371 ]. By contrast, if J,j < for some then 
the existence of solutions of Eq. © is dependent on the 
topology. 

The <&-notation allows to investigate which combina- 
tions of negative JV, in a graph Q lead to the violation of 
at least one of the Eqs. ©: The definition of $ implies 
that any given ^s^s can be represented as a function 



over a set of $ 
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where the Si arc subsets of S, and / depends on the 
choice of the Si. The key point is now that for show- 
ing that a given topological arrangement of negative Jij 
is incompatible with the stability conditions Eqs. ©J it 
suffices to show that there exists some ^si.s and some 
expansion (fT0|) thereof, for which there is no solution with 
all occurring <f> > 0. Let us emphasize that the sketched 
approach dispenses with the need to specifically calcu- 
lated any ^s^s, an d hence with addressing the up to 
l^l! terms they subsume. 

In [24[ , we employ the outlined approach to show that 
stability requires that every component of Q has a span- 
ning tree, whose links all have positive weights. A de- 
tailed mathematical proof of this result will be published 
in (2f| . For the sake of completeness we present a sketch 
of a proof in Appendix A. The proof uses the ^-notation 
to show that at least one of the necessary conditions must 
be violated when no spanning tree of positive connections 
exists. It thereby also provides an illustration of the di- 
rect application of the notation in a calculation. 

In addition to the spanning tree criterion for stability, 
which pertains to a global property of Q, Eq. ^ implies 
further restrictions, which pertain to mesoscalc proper- 
ties of Q. Thus, the weights of all links that are not 
part of any cycle of Q have to be positive. Moreover, at 
most one of the links that build an unbranched segment 
of a cycle of Q may have a negative weight. Finally, the 
weight of a 'negative link' in an unbranched segment of a 
cycle of G is bounded below by a value that depends on 
the weights of the other links in the segment (cf. Fig 2]). 

Note that while the mesoscale criteria restricting the 
number and position of negative links can be subsumed 
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under the global spanning tree criterion, the mesoscalc 
criteria restricting the weights of possible negative links 
arc inherently bound to the mesoscale. In particular, 
the latter can only be derived by considering ^s^s with 
small and intermediate \S\. This highlights the benefits 
of the proposed approach: The JSC provide stability cri- 
teria on all scales \S\, which are made accessible by means 
of the graphical notation and symbolic calculus. 

Applied to the Kuramoto model defined in Eq. ([TJ, 
the spanning tree criterion for stability allows to deduce 
properties shared by all possible phase- locked states [lij . 



In a phase-locked state, = Ay cos(xj 



where ; 



Xi is the stationary phase difference between oscillator j 
and i. Given that all link weights Ay > 0, stability of 
the phase-locked state thus requires that the coupling 
network has a spanning tree of oscillators obeying \xj — 
Xi\ < tt/2. 



V. GENERAL SYSTEMS OF COUPLED UNITS 

The derivation of stability criteria by graphical inter- 
pretation of the JSC is applicable to all systems with Hcr- 
mitian Jacobian matrix. Further, the simplification lead- 
ing up to the spanning tree criterion is possible, whenever 
the Jacobian has zero row sums. This condition is satis- 
fied for instance by all systems of the form 



Ci 



O l3 {x - Xl ), Vtel...JV, (II) 



where the Aj are the weights of a symmetric interaction 
network and the Oij odd functions. We emphasize that 
the approach remains applicable in heterogeneous net- 
works containing different link strengths Aj, coupling 
functions Oij or intrinsic parameters Ci. 

The class of systems to which the present results are 
directly applicable include general networks of phase os- 
cillators as well as other models such as continuous-time 
variants of the Deffuant model of opinion formation [4(| 
or a class of ecological metapopulation models |4lj . 

Although the specific criteria derived above are contin- 
gent on the zero-row-sum condition, it can be expected 
that the general approach proposed here is also applicable 
to situations where this condition is violated, such as in 
the model of cooperation among interacting agents stud- 
ied in [32| • A simple extension of Eq. (|TT|) which violates 
the zero-row-sum condition is found by replacing Ci by a 
function of Xi , which allows dynamical retuning of the in- 
trinsic frequency, e.g. for modeling homeostatic feedback 
in neural networks. For illustrating the application of the 
proposed method to models, in which the zero-row-sum 
condition is violated in some rows, we consider another 
model inspired by ncuroscience in the subsequent section. 



VI. ADAPTIVE KURAMOTO MODEL 

We now apply the proposed approach to an example of 
an adaptive network, i.e., a system in which the topology 
of the network coevolves with the dynamics of oscillators 
[26l [42|, 43] . In the context of the Kuramoto model, adap- 
tive coupling has recently attracted keen interest as it al- 
lowed to show that the emergence of synchronous motion 
can be intimately related to a selection mechanism of spe- 
cific network topologies (HI, and to the identification of 
complex hierarchical structures in the graph connectivity 

HE S3. 

We consider a system of N phase oscillators that evolve 
according to Eq. (QJ, while the coupling strength Aj 
evolves according to 

d 



dt 



Aj.j - — COSyXj — -^i) b ■ Ai 



(12) 



The first term in Eq. (|T2"T) states that the more similar 
the phases of two nodes the stronger reinforced is their 
connection, the second term guarantees convergence. In 
a stationary, phase-locked state state Aj = cos(xj —Xi)/b 
and all oscillators oscillate with a common frequency fi = 



JV Ei 



The stability of this state is governed by a 



symmetric Jacobian 
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where Oi 



j^i °ij i s ji 



\ cos 2 (xj - Xi), mi : 
sin(ccj — Xi) and we have chosen N = 3 for illustration. 
The marked partitioning separates two blocks on the di- 
agonal. The upper one is a diagonal submatrix of size 
L x L, L := N(N — 1)/2, the lower one is a N x N sym- 
metric submatrix with zero row sum, which we denote as 
j- 

Let us start our analysis by focusing on the upper left 
block of J. In the chosen ordering of variables, the first L 
minors D\s\,s satisfy the stability condition Eq. ^ iff b > 
0. Concerning the minors of order |5| > L, the following 
conventions prove advantageous: Below, we consider sets 
S that contain all variables Aj and n of the variables 
Xi. For every such set, we define S' as the subset of 
S that only contains the n variables Xi. Moreover, we 
define j as the matrix that is obtained from j if all are 
substituted by cos(2(a;j — Xi)). Finally, we define D n .s> as 
the determinant of the submatrix of j, which is spanned 
by the variables S S'. 

We find that 
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L+n,S 



(-l) L b L - n -D n „ 
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As j is symmetric and has a zero row sum, its minors, 
D n ^s, can be rewritten using Eq. (JT)) 
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where & nt s' refers to subgraphs of the graph Q defined 
by the off-diagonal entries of j. Stability requires that 
sgn(D^ +Ji s) = sgn (( — l) L+n ). As the necessary stabil- 
ity condition b > implies b L ~ n > 0, it follows that in a 
stable system 

>0, VS, n=l,...,rank(j). (16) 

Comparison with Eq. ((9]) reveals that a necessary condi- 
tion for stability is that every component of the graph Q 
has a positive spanning tree. Revisiting the definitions of 
graph Q, we find that the weight of a link ij of Q is given 
by cos(2(x :( — Xif). Hence, every component of Q has a 
positive spanning tree iff every component of the adap- 
tive coupling network has a spanning tree of oscillators 
obeying \xj — Xi\ < it/ A. The restriction on the station- 
ary phase-differences in a stable, phase-locked state are 
thus more strict in the adaptive than in the non-adaptive 
case. 



VII. APPLICABILITY OF THE GRAPHICAL 
NOTATION TO ISOSPECTRALITY PROBLEMS 

In addition to stability analysis, the graphical notation 
introduced here allows for instance exploring the isospec- 
trality of Hermitian or symmetric matrices p8l . |39| , which 
differ with respect to the signs of some off-diagonal en- 
tries. To see this consider a Hermitian matrix A G C nxn . 
Its characteristic polynomial x can be calculated as 

X (A) = A* (A - AI) 

Considering the structure of the graph Q associated to 
A — AI allows to determine which symbols contribute to 
X- One can then determine which changes of sign of off- 
diagonal entries leave all contributing symbols and thus 
the characteristic polynomial and the spectrum invariant. 

For instance, if the graph Q is a tree, then the only 
symbols that contribute to the characteristic polynomial 
X are the symbols x and |. Thus, % is a polynomial 
of factors (An — A) (corresponding to symbols x), and 
factors AijAji (corresponding to symbols |). Due to the 
hcrmiticity of A, A^Aji — |Ay| 2 . It follows that Xi and 
therewith the spectrum, is invariant under any operation 
that changes the sign of a pair of off-diagonal entries 

A'ij y ^iji ^-ji ^ ^-ji' 

Along the same line, we can infer isospectrality rela- 
tions for matrices A, whose corresponding graphs Q are 
composed only of tree-like subgraphs and isolated cycles. 
For such matrices, the spectrum is invariant under sym- 
metry preserving sign changes of 

• an arbitrary number of off-diagonal entries that do 
not belong to cyclic subgraphs of Q. 

• an even number of off-diagonal entries that belong 
to the same cyclic subgraph of Q. 



VIII. CONCLUSIONS 

In the present paper we presented a general approach 
for the analysis of symmetrically coupled systems and a 
specific simplification for systems, in which the Jacobian 
matrix J additionally has a zero row sum. 

Using a graphical interpretation of Jacobi's signature 
criterion we showed that for all non-trivial eigenvalues of 
J to be negative, the graph Q, whose adjacency matrix is 
given by the off-diagonal part of J, has to obey necessary 
topological conditions. Thus, Q must have a spanning 
tree of links with positive weights, which restricts the 
number and position of potential negative entries in the 
Jacobian matrix. Moreover, the absolute value of poten- 
tial negative link weights is bounded above by a topology 
depended relation. 

We used the approach for deriving necessary conditions 
for local asymptotic stability of stationary and phase- 
locked states in networks of phase oscillators. Our results 
provide an analytical angle that is complementary to sta- 
tistical analysis of network synchronizability. Where sta- 
tistical approaches reveal global features impinging on 
the propensity to synchronize, our approach can pinpoint 
specific defects precluding synchronization. We note that 
such defects can occur on all scales, corresponding to the 
violation of the signature criterion in subgraphs of differ- 
ent size. This highlights synchronization of phase oscilla- 
tors as a simple but intriguing example in which instabil- 
ities can arise from local, global or mesoscale structures. 
In the future, the approach proposed here may provide a 
basis for further investigation of these instabilities. 

A limitation of the present approach is that the Ja- 
cobian matrix of the state under consideration has to 
be available in some form. However, this limitation is 
inherent to all forms of linear stability and bifurcation 
analysis and several strategies have been developed for 
mitigating it. These include for instance the equation- 
free approach (45j . where the desired information on the 
Jacobian is generated on-demand from simulations or ex- 
periments, and generalized modeling |46| . which reveals a 
parameterized presentation of the Jacobian for very gen- 
eral classes of systems. Even if these specific techniques 
are not applied, researchers are often aware of the struc- 
ture of the Jacobian in the system under consideration, 
which can be sufficient for gaining some fundamental in- 
sights by the proposed approach. 

It is apparent that applying the stability analysis pro- 
posed here may require overcoming specific problems. In 
particular, the requirements of symmetry and zero row 
sum may necessitate extending the basic scheme or focus- 
ing on specific cases, as we did in the adaptive network 
example. However, if these problems can be overcome, 
the approach can reveal hard analytical conditions that 
all stable steady states must meet. Further, it can iden- 
tify topological structures (on all scales) that prevent the 
existence of such states and pinpoint them in networks. 
Because this type of information is complementary to 
insights revealed by common statistical approaches, it 
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potentially of high value for the analysis of the system. 



Appendix A: Spanning Tree Criterion for stability 

To prove the spanning tree criterion for stability con- 
sider that in a network without a positive spanning 
tree it must be possible to partition the nodes into two 
nonempty sets 1%, I2 such that 



Jij < V i, j I i e h, j e h 



(Al) 



The idea is now to evaluate the stability conditions © 
for different S C I 1 thereby exploiting that all links lead- 
ing out of I\ have negative weights. For this purpose, it 
is convenient to define E* as the set of links connecting 
Ii and I 2 , and X = {x\, . . . ,x m } as the subset of nodes 
€ I\ incident to at least one link from E* ('boundary 
vertices'). Further, we define Ui as the sum over all ele- 
ments of E* incident to x i; and, for any subset Y of X, 



cry 



n 



rnCY 



a m . Finally, for any subset F of 1, we 
define ry as the sum over all forests of Q that (i) span 
I\, and (ii) consist of \Y\ trees each of which contains 
exactly one element from Y. 

With the above definitions we can write 



BCX\C 



(A2) 



where B and C are disjoint subsets of X and 



q,h\C- 



The rather abstract equation is illustrated in 



Fig. [5] by means of an example. 

Equation (|A2[) can be read as an expansion of ^i ± \c 
in contributions from E* . This expansion has the advan- 
tage that the sign of one of the two factors in each term 
is known: sgn(oB) = (— l)' s '- The signs of the factors 
t buc- however, are still indeterminate. Nevertheless, we 
can show that Eq. (|A2[) is incompatible with the stabil- 
ity condition ([9]) by considering a linear combination of 
ii\C- 

M) |C V $/AC = (- 1 ) |C|(J C ABTBUC 

BCX\C 



ccx 



>0 by >0 
construction 



CCX 



\\C\ 



CCX 
BCX\C 



ACX 
CCA 



ACX CCA 



which is a contradiction. Therefore the existence of a 
spanning tree of positive elements is a necessary condition 
for stability. 

Appendix B: Lower bound on the negative link 
weights 

Consider a path of d — 1 degree- two nodes t)j, i 6 
{1, . . . , d — 1}, that are part of at least one cycle of Q. 
Together with the d egdes Ci that are incident to at least 
one of the nodes Wj, the path constitutes an unbranched 
segment C of a cycle of Q (cf. Fig. 3). 

According to the spanning tree criterion, C can have at 
most one link with negative weight. Below, we consider 
the case where C has exactly one such link, and show 
that Eq. (|9]) imposes an upper bound on the absolute 
value of the negative link weight. 

We use the following conventions: Let the nodes be 
labeled such that the indices i occur in an increasing 
order if C is paced out. And let the links be labeled such 
that i>i is incident to c, and q+i (cf. Fig. 2S). Further, 
let wi denote the weight of c;. And lastly, let c x denote 
the only link in C with negative weight w x . Below, we 
show that stability requires that 

\ w I < flail Wi, i g J* 

I x I ) ' all distinct products of (d — 2) factors Wj, i g J* ' 

(Bl) 

where I* — {1, . . . , d} \ x. 

For deriving Eq. (|B1|) . we consider a sequence of 
conditions <£% . . . &s d _ 1 ■ The sequence is constructed 
as follows: we choose Si = x, S2 = {x,x + l}, 
S3 = {x, x + 1, x + 2}, and so forth until Sd- X = 
{x, x + 1, x + 2, . . . , d — 1 }. The remaining elements of 
the sequence are then constructed as S^-x+i = Sd- X U 
{x - 1}, S d ~x+2 = Sd-x U {x— l,x — 2}, and so forth 
until S d -i = {«i} ie{ i,..., d -i}- 

The first element of the sequence, $^ > 0, stipulates 
that w x + w x +i > 0, and thus that —w x < w x +i- 

The second element of the sequence > stipulates 
that WxW x +i + w x w x +2 + ui x +\w x +2 > and thus that 

-W x < (Wx+lWx+2) / (w x +l +W1+2). 

More generally, every element of the sequence 
$5 15 . . . , $s d _ 1 modifies the upper bound on — w x as per 



-w x < 



Y[ all Wi, i £ S* 
all distinct products of (\Si \ — 1) factors wi, i G S i 



=0 



(B2) 

where S* = (Si \x)U (max(5f) + 1). 

The right hand side of Eq. (|B2I) is monotonously de- 
creasing with increasing \Si\. We can thus conclude that 
C can have at most one link with negative weight, whose 
absolute value \w x \ is bounded above by Eq. (|B 1[) . 
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